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PROCEDE DE MESXJRE DE PARAMETRES PHYSIQUES D'AU MOINS 
UNE COUCHE D'UN MATER1AU AUX DIMENSIONS MICROMETRIQUES 

La pr^sente invention a pour objet un procMi de determination d'au 
moins une caractdristique mdcanique d*au moins une couche disposde k la surface 
5 d*un substrat, notamment cristallin. 

La connaissance des propriety m^caniques d*un syst^e ou 
dispositif physique permet d'en optimiser le fonctionnement. Tout mat^au quel qu*il 
soit est soumis k des contraintes extemes et il faut pouvoir conna!tre sa resistance a de 
telles contraintes. 11 est done important de connaitre au mieux les propri6t6s elasliques 
10 d*un sysl^me et en particulier cellos d'une ou plusieurs couches disposes h la 
surface d'un substrat. Les propriety de cette couche (ou de ces couches consider^es 
comme un sous-systeme) different en effet notablement des propriSt^s des couches 
prises individuellement. 

Les propri6t6s 61astiques des maieriaux interviennent dans de 
15 nombreux domaines d'applications : revetement de pieces m^caniques, deformation de 
structure, ... EUcs jouent un role de plus en plus important dans la realisation de 
circuits electroniques : la taille des dispositifs diminuant, les contraintes generdes k 
Tinterface des difT^rentes parties du dispositif sont d'autant plus importantes. 

Les methodes d^ja developp^es pour mesurer les constantes 
20 d'eiasticite et les coefficients de dilatation sont g6neralement realisdes a une echelle . 
macroscopique. 

Les methodes traditionnelles sont : 

A) deflexion, deformation d*eprouvettes macroscopiques ou de taille 
plus reduite : voir par exemple : Measurement of elastic modulus, Poisson ratio, and 
25 coefficient of thermal expansion of on-wafer submicron films. Zhao, Jie-Hua; Ryan, 
Todd; Ho, Paul S.; McKerrow, Andrew J.; Shih, Wei-Yan. PRC/MER2.206, 
Journal of Applied Physics (1999), 85(9), 6421-6424, 

ou bien: Elastic modulus measurements by three methods on a 
2-D laminated continuous fiber-reinforced ceramic matrix composite. Gonczy, S. T.; 
30 Jenkins, M. G Ceramic Transactions (2000), 103(Advances in Ceramic-Matrix 
Composites V), 541-547, 

ou bien encore: Analysis of residual stress in cubic boron nitride 
thin films using micromachined cantilever beams. Cardinale, G. F.; Howitt, D. G.; 
McCarty, K. F.; Medlin, D. L.; Mirkarimi, P. B.; Moody, N. R. Diamond and Related 
35 Materials (1996), 5(11), 1295-1302 



2 



B) Ondcs 61ectroacoustiques : par cxemple Hardness and Young's 
modulus of high-qualily cubic boron nitride films grown by chemical vapor 
deposition. Jiang, X.; Philip, J.; Zliang, W. J.; Hess, P.; Mntsumoto, S Journal of 
Applied Physics (2003), 93(3), 1515-1519, ou bien Thin-film elastic-property 

5 measurements with laser-ultrasonic SAW spectrometry. Hurley, D. C. ; Tewary, V. 
K.; Richards, A- J Thin Solid Films (2001), 398-399 326-330 

C) nanoindentation par exemple : Determination of mechanical film 
properties of a bilayer system due to elastic indentation measurements with a spherical 
indenter. 

10 Chudoba, T.; Schwarzer, N.; Richter, F.; Beck, U. Thin Solid Films 

(2000), 377-378 366-372, ou bien. Hardness and elastic modulus measurements of 
AIN and TIN sub-micron thin films using the continuous stiffness measurement 
technique with FEM analysis. Rawdanowicz, T. A.; Sankar, J.; Narayan, J.; 
Godbole, V. Materials Research Society Symposium Proceedings (2000), 594(Thjn 
1 5 Films-Stresses and Mechanical Properties Vlll), 507-51 2. 

11 existe 6galement des m6thodes plus r6centes pour mesurer les 
propri6t6s m6caniques de petits objets , par exemple la microscopic k force atomique 
(•AFM*) d6crite par exemple dans : 

Measurements of elastic properties of ultra-thin diamond-like carbon 
20 coatings using atomic force acoustic microscopy. Amelio, S.; Goldade, A. V,; Rabe, 
U.; Scherer, v.; Bhushan, B.; Arnold, W Thin Solid Films (2001), 392(1), 75-84. 

L'invention a pour but de mesurer un ou plusieurs parametres 
m6caniques, notamment les constantes £lastiques et les coefficients de dilatation 
thermique d'un materiau et les contraintes que ce materiau g6nere lorsqu'il est associe 
25 a d'autres. Sa sp6cificit6 est de pouvoir mesurer de telles constantes dans des 
conditions proches de leur condition d'applications : couches minces ou dpaisses, 
couches inhomogenes, 'couches' discontinues ou pr6cipites/inclusions, boites (terme 
plus specifiques aux matdriaux semi-conducteurs ou les charges sont local isees dans 
ces boTtes constituees d'un second materiau different du substrat). Les connaissances 
30 actuelles ne permettent pas d'affirmer que les constantes mesurdes 
macroscopiquement peuvent fitre utilis^es dans des structures de petites dimensions, 
c'est a dire micrometriques et nanom6triques. 

Uid6e de base de Tinvention consiste a utiliser un ph6nomene bien 
connu, mais qui est g6n6ralement considere comme un inconvenient que Ton cherche 
35 a 6viter ou a negliger : dans une lame mince, du fait de la proximit6 de la surface libre. 
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les contrainles accumul^es dans le dispositif initial sont relaxes (voir les documents 
pr6citds dc Perovic, D., Weatherly G. and D, Houghton Phil. Mag. 64 (1991)). 

Le dispositif & ^tudier doit comporter deux parties (figures la & Id) : 

-a) un substrat A de r6f<§rence k propri^tds m^caniques bien connues 
5 (par exempic du silicium, du germanium, du saphir» un syst^me m^tallique ...) 

-b) une ou plusieurs parties B & la surface de ce substrat de rdf^rence. 

La prdsente invention consiste k : 

(i) pr61ever une lame mince k g6om<itrie contr616e dans le sysl6me 
£tudi6. L'6paisseur de cette lame mince doit en particulier 6tre suffisamraent faible 

10 pour que les contraintes au soin de cette lame soient effeclivement relSch^es. 

(ii) mesurer les deformations de cette lame. Diff6rentes 
temperatures doivent etre utilis6es si Ton veut determiner precis^ment les coefficients 
de dilation thermique. 

(iii) dventuellement r6amincir la lame pour r6p6ter les actions (i) et 
1 5 (ii) avec une lame d'6paisseur differente. 

(iv) utiliser les valeurs mesur6es des deformations de ces lames de 
differentes ^paisseurs pour determiner en particulier les constantes eiastiques et lesi. • 
coefficients de dilation des differents materiaux de la partie B, notamment en utilisant . 
une technique de modeiisation. 

20 Selon Tinvention, il tfest pas n6cessaire de faire des mesures directes;- 

sur le materiau a caracteriser (partie B) pour determiner les constantes de ce materiau: 
II reste bien sur possible de r^aliser des mesures egalement sur la partie B et cela 
apporte un complement d'information. Ce sont les deformations/rotations du substrat 
(partie A), mesurees quantitativement qui pemiettent de remonter aux caracteristiques 

25 des parties B. Ainsi, on determine les proprietes de la partie B en analysant la fa9on 
dont la partie A se deforme dans une lame amincie ou dans des lames de differentes 
epaisseurs. 

Dans le systeme le plus simple, la partie B est composee d'une seule 
couche mince ou epaisse deposee sur le substrat. Des systemes plus complexes 
30 (plusieurs materiaux, gravures, petites particules) peuvent egalement etre analyses, 
mais ils necessitent un plus grand nombre de parametres a minimises 

L'invention conceme ainsi un procede de determination d'au moins 
une caracteristique mecanique d'au moins une couche disposee k la surface d^un 
substrat caracterise en ce quMl comporte : 
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a) la realisation d*une lame L d'ipaisseur t suffisamment faible et 
ayant deux faces sensiblement paralldles et disposics sensiblement 
perpendiculairement h ladite surface de substrat ; 

b) la mesure sur ladite lame d'au moins un param&tre de 
5 ddformation du substrat h difPferentes profondeurs par rapport a la surface; par 

parametre de deformation, on entend ^galement paramitre de rotation du substrat ; 

c) la determination h partir au moins dudit parametre de 
deformation/rotation d'au moins une caract6ristique m6canique de ladite couchc. 

Le procide peut comporter la realisation de plusieurs lames 
10 d^epaisseurs differentes ainsi que la mise en osuvre de retape sur chacune desdites 
lames. Pour au moins une dite lame, retape peut 6tre repetee h au moins deux 
temperatures differentes. 

Ladite mesure est avantageusement realisee en generant, pour, des 
points du substrat situes a differentes profondeurs, des diagrammes de diffraction d'un 
1 5 faisceau eiectronique convergent d'axe Zo desoriente par rapport a la normale a ladite 
lame, lesdits diagrammes comportant des lignes ou bandes de Holz. La determination 
c peut alors comporter le releve de la largeur des lignes de Holz d'au moins certains 
desdits diagrammes, pour au moins un plan cristallographique du substrat. A partir de 
la largeur de ces lignes de Holz, on peut calculer pour chaque diagramme une rotation 
20 maximale Betamnx le long de Taxe du faisceau d'electrons. Cette rotation est induite 
par la couche (ou les couches) disposee sur le substrat et elle en caracterise les 
proprietes. 

On peut ensuite tracer au moins une courbe representant ladite 
rotation Beta^ax en fonction de la profondeur a laquelle ont ete obtenus lesdits 

25 diagrammes.On peut ensuite tracer par simulation au moins une courbe repr6scntant 
cette rotation en fonction de la profondeur pour des valeurs possibles des coefficients 
a extraire par exemple le module d'Young et/ou le coefficient de Poisson de la couche 
a analyser dans le cas d'une approximation isotrope et minimiser recart entre au 
moins une courbe simuiee el une courbe experimentale correspondante, pour 

30 determiner ces coefficients. Une technique similaire pourra etre utilisee dans le cas 
d'unemodeiisation anisotrope avec les coefficients connus de THomme du Metier. On 
prendra pour valeur les parametres ayant conduit a la courbe simuiee la plus proche de 
la courbe experimental^ 

LMnvCTition sera mieux comprise k la lecture de la description ci- 

35 apres, en liaison avec les dessins dans lesquels : 
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- les figures la & Id illustrent difT<§rentes configurations du dispositif 
k dtudier : couche simple pour la figure 1 a, bandes de metallisations ou Hots pour la 
figure lb, couclie avec zone incluse dans le substrat pour la figure Ic, et transistor 
pour la figure Id. 

5 - La figure 2a et 2b repr^entent une vue en coupe d'un transistor 

ayant un subslrat en silicium Si, surmont6 d'une couche de NjSi d'^paisseur 
e = 20 nm et d'une couche de Si3N4. Un contact en tungst^ne constitue le contact D de 
drain. La zone analys6e (figure 2b) est point6e par la fleche F. La lame mince 
representee a la figure 2b prdsente une 6paisseur t (que Ton peut faire varier). Le 

10 faisceau incident moyen des electrons pour les cliches "CBED" (direction zo) est pris 
selon Taxe ya qui fait un angle y avec Taxe yi normal k la lame. Les directions et les 
axes cristallographiques correspondant sent illustres a la figure 2b- 

La figure 3a est un montage de photos illustrant cinq des 
digrammes de diffractions "CBED*' choisis panni une cinquantaine dMmages realisees 

IS effectivOTient, selon unedroite perpendiculaire a la surface (direction - Z2)» a une 
distance de 1 55 nm du contact D de drain. Pour chaque diagramme, on indique la 
distance - Z2 (en nanometres) par rapport a la surface du substrat (interface avec la ^ 
couche mince) et Tangle AO mesur6, la courbe de la figure 3b representant revolution .... 
de Tangle A6 calcuie en fonction de -Z2. ^ 

20 - La figure 4a montre un cliche de diffraction « CBED » pris . 

suivant un axe de zone [230]dans le silicium non contraint, a une profondeur 
importante Z2 = 300 nm (X2 // [3, -2, 0], y2, // [2, 3, 0], Z2 //[OOl ]. 

La figure 4b est le cliche 4a sur lequel on a ajust6 des lignes de 
Holz positionnees par le logiciel "jems". La position de ces lignes permet de 

25 determiner avec precision Tangle de demi-convergence du faisceau incident et Ja 
tension d*acceieration (la methode est ici similaire a celle utilisee dans le projet 
"STREAIVT). 

La figure 5a est un cliche CBED en axe de zone [230] dans le 
silicium defoime (diagramme « CBED » de type c pris k la profondeur Z2 = 1 40 nm). 
30 - Sur la figure 5b, on a superpose a ce cliche des lignes issues de 

3 simulations avec le logiciel " jems " : 

■ le prCTfiier systeme de lignes rq[)r6sente la simulation du cliche 
de diffraction d'un cristal de silicium parfait desoriente de +Omflx par rapport a Taxe 
X2// [3-20] ; 
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- le deuxi^me syst&me de lignes repr6senle la simulation du 
clich6 dc diffraction d*un crista! dc silicium parfait d6sorient6 de -Gmax par rapport h 
Vaxc X2//[3-20]; 

le trojsieme systfeme de lignes (en poinlilltes) rcpr6sente la 
5 simulation du clich6 de diffraction du cristal parfait non d6sorient6. On note par AO, 
Pangle 2G,nax.. 

- Les figures 6a k 6d illustrent la faisab)lit6 de minimiser les 
coefficients d*elaslicit6 du matdriau. 3 parametres ont 6t6 minimises manuellement par 
essais/erreurs : la temperature de cohdrencc To (temperature fictive ou les deux 

1 0 materiaux seraient coh6rents), le module d^Young E et les coefficients de Poisson v. 

Les cairns repr6sentent les mesures exp6rimentales de Tangle A9. Les 
courbes avec de nombreux points sont obtenues a partir des simulations par elements 
finis et repr6sentent Tangle 0,98AP ou AG = 2 Beta^^n,,. Le critere de minimisation % 
6value la distance entre les angles AG mesur6s et 0,98AP calculus, a diflKrentes 
1 5 profondeurs zo sous le siliciure et h diff^rentes 6paisseurs t de lame. 

La figure 6a : 6paisseur de la lame t = 300 nm, 
valeurs initiales des parametres pris dans la litt6rature : T = 410*'C, 
E= 150 GPa,n = 0,l 

La figure 6b : 6paisseur de la lame t = 300 nm, 
20 valeurs initiales des parametres pris dans la litt6rature : To = 410**C» 

E= 150 GPa,n = 0,l. 

La figure 6c : dpaisseur de la lame t = 320 nm, 
valeurs finales apres une minimisation manuelle 
partielle : To = 430^C, £-115 GPa, n = 0,288. 
25 La figure 6d : 6paisseur de la lame t = 320 nm, 

valeurs finales apres une minimisation manuelle 
partielle : To = 430^C, E = 115 GPa, n = 0,288. 

La figure 7a est un diagramme de diflfraction CBED obtenu 
pour une lame d'epaisseur t = 320 nm a la distance -22 = 84 nm de la surface. 
30 - La figure 7b est un profil qui a 6t6 realist sur la bande de Holz 

(5,-3,7) de la figure 7a. L'6cart angulaire AGg est marqu6 par 2 traits verticaux et 
repr6sente la largeur de la bande de Holz. 

La figure 7c est une simulation du profil de la figure 7b qui 
montre quMl est possible de reproduire T^largissement des lignes de Holz et les 
35 variations d'intensitd 1 - (9) dans les bandes de Holz. Cette simulation a utilise les 
r6sultats du calcul par Elements finis qui est une illustration du point iv (constantes 
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dastiques optimises et d^placement R(ywZ|)). AOg est 6gal a 0,1 1"^, ce qui correspond h 
une rotation d'axe X2 d*angle A6 - 0,14^. Dans les calculs prdliminaires, on n'a pris en 
compte que les largeurs ABg d'axe X2//[320]. Un algorilhme de minimisation plus 
perfectionni chercherait h reproduire les courbes similaires h 7b dans leur ensemble, 
avec leur oscillation et non pas h reproduire seulement leur largeur A9g h travers 
Tangle AG. 

Les figures 8 & 8d illustrent un deuxiime exemple experimental, 
dans lequel le substrat A de silicium est surmont6 d'une couche de Si(i.x) Gex, puis 
d*une couche de Si en surface. La figure 8a montre la lame fine d*^ai$seur t. Les 
figures 8b a 8d sont des diagrammes de diffraction " CBED " respectivement dans la 
direction [230] a I'intdrieur de la couche Si(i.x) Ge^, dans le substrat Si loin de la zone 
d^form^e, et dans le substrat Si prfes de la zone d6form6e. 

On pr6sentera ci-aprte les quatre points pr6cit6s de IMnvention. 

point (i) : La lame mince h g6om6trie contr616e est extraite ou 
amincie dans le dispositif. Une lame a face paralldle est pr6f§rable mais pas 
indispensable. Un 16ger angle peut 6tre present. On a utilise un faisceau d'ions focalis6 
TIB', mais des procedes altematifs, classiques en preparation d'echantillons pour la 
microscopic dectronique peuvent etre utilis6s (amincissemeni mdcanique, clivage,...) 
mais la technique TIB' a Tavantage d'€tre rapide et de ne pas perturber 
m6caniquement le systSme ou dispositif et de contr61er entiferement les operations. 

point : Pour mesurer les deformations, on utilise par exemple un 
faisceau d'eiectrons convergent (CBED) en mode balayage. Un faisceau eiectronique 
de taille nanometrique est focalise en differents points de rechantillon de fa9on a 
cartographier les deformations completes de la lame mince (figure 3). En. chaque point 
un diagramme de diffraction 'CBED* est done obtenu. Ce diagramme comporte de 
nombreuses lignes de Holz (plus de 1 0 lignes) qui correspondent chacune S un plan 
cristallographique indexe par un vecteur g du reseau reciproque. 11 est preferable 
d'acquerir les diagrammes 'CBED' dans le substrat (partie a), car la partie b est 
general ement trop fine ou trop defectueuse pour que Ton puisse acquerir de bons 
cliches de diffraction dans cette zone b. Mais pour certains systemes bien cristallises, 
des cliches peuvent etre obtenus dans la zone a et b. Ceci serait par exemple le cas si 
la couche a est du silicium et la couche b est composee de deux couches fines, une 
couche de Sij.xGex (composition x faible en Ge, par exemple 10%) et une couche 
de Si. 



Pour une 6paisseur de lame donn^e et une temperature d'observation 
donn^e tout un ensemble de cliches CBED est rdalisi de fa9on h d^temiincr les 
rotations et param&tres cristallins du substrat. 

point fiii^ : En faisant varier T^paisseur de lame (la technique 'FIB* 
5 pcrmet de r6aliser cela facilement) et la temperature d'observation dans le microscope, 
on oblient tout un ensemble de donn^es exp6rimentales qui permettront (point iv) de 
calculer les constantes d'dlasticitd, le coefficient de dilatation de la partie B ainsi que 
les contraintes dans 1c dispositif initial non aminci. Dans le cas ou la partie b est 
compos6e d\ine couche homogene d'un mat^riau donn6, il n'est pas nteessaire que les 
10 mesures soient r6alis6es sur une m€me lame successivemenl amincie h difftrentes 
dpaisseurs. Travailler sur une lame unique augiTiente toutefois la precision et est 
indispensable dans le cas ou le systfeme est constitu6 d'un nanosyst^me (transistor par 
exemple) unique. Une seule 6paisseur de lame ne pemiet pas de calculer correctement 
loutes les constantes du materiau. Une seule 6paisseur de lame nc donne des 
15 informations que sur les contraintes du mat6riau (c'est le cas de I'^tude de la courbure 
des tranches semiconductrices ou ^wafers' via la fonnule de Stoney). Les contraintes 
sont en partie rel&chees par une courbure du substrat (voir par exemple : Measurement 
of elastic modulus, Poisson ratio, and coefficient of themial expansion of on-wafer 
submicron films. Zhao, Jie-Hua; Ryan, Todd; Ho, Paul S.; McKenrow, Andrew J.; 
20 Shih,Wei-Yan. Journal of Applied Physics (1999), 85(9), 6421-6424) 

point (iv't : Des simulations utilisant la th^orie ^lastique sont ensuite 
effectuees pour reproduire les rteultats exp^rimentaux. Dans des systftmes complexes 
seuls des calculs par 616ments finis peuvent 6tre r6alis6s. Dans des systdmes plus 
simples (partie B = couche mince) une formule analytique peut etre impl6ment6e. Les 
25 simulations reproduisent les ph6nomenes suivants : 

- Dans le dispositif initial, les differenles parties sont contraintes ou 
partiellement contraintes - Le fait d'extraire une lame mince du dispositif initial 
relache les contraintes sous la fomie de rotation/changement de parametres cristallins 
et c'est Tobservation et la simulation de cette relaxation des contraintes qui nous 
30 permet de determiner les caractdristiques de la partie B. 

L'invention trouvera de nombreuses applications dans le traitement 
de surface de pieces m6caniques, d'optimisation de circuits electroniques (contact 
mdtallique, couche d'oxyde ...) ou dispositifs ou la presence de deux materiaux 
differents cr6e necessairement des contraintes mecaniques. 
35 Le proc6d6 selon I'invention est original bien qu'il fasse appel a des 

techniques ou des effets physiques bien connus : 
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la technique du faisceau d'ions focalis^ ou TIB* (Focus Ion 
Beam) pour la pr^aration d'ichantillon. 

le faisceau convergent (CBED : Convergent Beam Electron 
Diffraction), une technique particuliere k la microscopic dectronique. Le projel 
5 STREAM (voir par exemple la publication 'Software for automation of TEM/CBED 
Methodology for strain determination' lST-1 9999-1 0341 STREAM Consortium - 
Deliverable D23) utilise une technique proche (TIB*, faisceau convergent, 
simulation), mais la technique mesure essentiellement des variations de param^tres 
cristallins alors que le proc6d6 selon I'invention s'int6resse principalement aux 
10 rotations locales du r6seau cristallin. De plus le projet STREAM ne cherche pas k 
mesurer des constantes 61astiques, mais k mesurer des deformations contraintes dans 
des circuits int6gres. 

Le projet STREAM ne mesurait que des variations de parametres 
cristallins loin des deux parties A et B du dispositif et n^gligeait la relaxation des 
15 contraintes dans la lame mince. La detection de rotation du r6seau cristallin selon 
rinvention permet d'etre plus rapide, plus precis et de se rapprocher de I'interface 
entre les parties A et B. 

La puissance du proc^de selon Tinvention a 6t6 montr^e en analysjoit 
les d6fonnations introduites par une couche de NiSi dans un circuit int6gr6 
20 (6chantillon n**l , figure 2a). Ce systfeme est relativement complexe, car il y a plusieurs 
contacts dlectriques, plusieurs mat6riaux pr&ents. Pour 6tie trait6 rigoureusement 
(ce qui est tout a fait possible), il faudrait tenir compte de toutes les composantes du 
systfeme. 

Dans une premiere analyse, on a suppos6 que la partie B.pouvait 6tre 
25 considdr^e comme une couche mince de taille lat^rale infmie surmonlte par 
Tatmosphere (les couches situ6es au dessus de la couche NiSi dtaient 
experimentalement decoll6es, I'^aisseur de la couche mince a 6t6 mesur6e 6gale h 
e=20nm (voir figure 2b). 

On d6crira ci-apres de fagon precise la realisation des quatre points 

30 de la partie 3. 

Pour cela, il faut ddfinir diff6rents rep^res g6om6triques. 

Ro =(xo,yo,Zo) ddsigne le repfere geom6trique lie au microscope. L'axe 
2o est paralieie h l'axe optique du microscope, defini comme etant la direction 
moyenne selon laquelle les electrons se propagent avant I'echantillon. Les images ou 
35 les cliches de diflraction sont enregistres dans le plan (xo,yo). 
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On d^finit dgalemenl difftrcnts repferes Rc'='(Xc,yc,2cX Ri=(xi,yi,Zi), 
R2='(X2,y2,Z2), li6s a la structure crislalline du subslrat dc silicium (partie A) : 

Xc=llOO), yc^lOlO], 2c=l001] 
5 X, =l/V^[lTO], y, =l/^/2(no], z,=(0011 

X, =1/Vi3[320], yj =l/>/l3[230], 22=[001] 

Le rcpere R2 se ddduit du repfere Ri par une rotation d'axe [001], 
d'angle T^l 1 .3 1*". (voir sch6n)a 2b). Les relations entre le repfere du microscope et les 
1 0 repdrcs cristallins dependent de rorientation de T^chantillon dans le microscope. 

Point (i^ : La lame mince a faces pratiquement paralleles a et6 
r6alis6e avec un TIB' (voir remarques 1 el 2). La normale aux faces a dt6 choisie trfes 
proche de la direction yi, (figure 2b) (mais ici aussi une gdometrie difn&rcnie pourrait 
€tre choisie : par exemple y2 ). Dans la premiere s^rie d'experience, T^paisseur t de la 
15 lame a 6le de t=320nm. Cette 6paisseur a 6te mesur6e par une technique relativement 
classique de faisceau convergent - voir le livre 'Electron MicroDiffraction* de J.CH. 
Spence et J.M Zuo (Plenum Publishing Corporation). 

Remarques : 1) L'6cart de parall61isme a 6t6 mesur6 par faisceau 
convergent et perte d'6nergie. Ces techniques sont relativement classiques (cf. le livre 
20 pr6cit6 'Electron MicroDiffraction' de J.CH. Spence et J.M Zuo pour le faisceau 
convergent et le livre 'Energy-Filtering Transmisssion Electron Microscopy de 
Reimer public chez Springer Veriag pour les pertes d'6nergie). On a done mesxirS 
I'^paisseur de la lame a difTArentes distances de la couche sup6rieure. On a tTOuv6 que 
la lame sup6rieure fait un angle de MS*" par rapport h la lame inferieure. Dans notre 
25 demonstration de faisabilit6, cet angle a 6x6 n6glig6 et on a considere que la lame avait 
des faces paralleles. Mais il aurait eX6 possible d'introduire cet angle dans les calculs 
d'616ments finis. 

2) Le FIB est 6galement capable de produire des lames avec un 
meilleur parallelisme. 

30 Point fii> : Les diagrammes 'CBED' ont ete pris dans (ou tres proche 

de) la direction d'orientation y2=[230], c'est a dire que y2 est parallfele a Zq. Cest la 
direction utilisde dans le projet STREAM, (mais d'autres directions d'observation sont 
egalement possibles, Uobtention dc diagrammes CBED dans diff6rentes directions 
viendrait augmenter le nombre de donn^es experimentales). Typiquement, la taille du 

35 faisceau 61ectronique a 6t6 prise egale a 0,4 nm. Tangle d'ouverture du faisceau est 
voisine de 15mrad. Lors d'un balayage, de 50 a 100 diagrammes CBED sont pris tous 
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les 4 nm dans une direction perpendiculaire h la surface (direction Z2), cn partant de la 
surface, mais seulement une dizaine est retcnuc dans les calculs. La figure 3 tndique la 
position scion Taxe Z2 ou S de ces cliches exp^rimentaux ont i{€ obtenus. 

Plusieurs types de diagrammes CBED peuvent Stre ddfinis : 
5 a) Loin de I'interface entre les parties A et B, le substrat de silicium 

est, avec la precision de mesure, consid£r£ comme non d6form6 (parexemple le 
diagramme -Z2='300nm de la figure 3a). 11 peut servir de r6f6rence pour determiner les 
param^tres exp^rimentaux comme cela est fait dans le projet STREAM. 

b) En se rapprochant de Tinterface, les lignes de Holz restent fines 
10 mais commencent h se d6placer 16gferement : les paramdtres cristallins du substrat 
subissent de tr&s I^gdres modifications. C'est Teffet mesur^ et quantifi^ dans le projet 
STREAM. 

c) Plus pres de Tinterface, les lignes de Holz s'dlargissent (par 
exemple le diagramme -Z2=139nm de la figure 3a). Nous les appelerons bandes de 
1 5 Holz. Cest Teffel qui est mesur6 et quantifi6 ci-apr6s. 

d) Encore plus proche de Tinterface, les lignes de Holz deviennent 
trop larges et trop faibles : les diagrammes *CBED* ne contiennent plus assez de 
details pour etre quantifies (par exempje le diagramme -Z2=70nm de la figure 3a). 
Point (iii^ : 

20 Une fois cette s6rie de mesures r6alis6e a I'epaisseur t=320nm, la 

lame a &t€ amincie une nouvelle fois jusqu'a une dpaisseur de t=^300nm. Une deuxidme 
s6rie de mesures a alors 6i& r^alisee. 

On a choisi 14 diagrammes pour la premiere Spaisseur et 10 
diagrammes pour la deuxieme (voir figure 5). Ce nombre est suffisant pour reproduire 
25 les variations de la rotation. 

Poinliiy): 

On a utilise le logiciel "jems" (P. Stadelmann, CIME-EPFL 
CHI 01 5- LAUSANNE) pour reproduire la position des lignes de H0J2 des cliches 
CBED loin de Tinterface (cliches CBED de type a), mais d'autres logiciels ou la 

30 th^orie presentee dans le livre pr6cite de Spence et Zuo pourraient etre utilises. Si Ton 
reporte ces lignes de Holz du substrat de silicium - le cristal de r6ference - sur les 
cliches CBED de type c, ces lignes se situent au milieu des bandes de Holz (la lame 
est homogene et symetrique). Les positions d'une des extremites de ces bandes de 
Holz sont simul6es dans le logiciel "jems" en ddsorientant le cristal de r^fdrence par 

35 une rotation d'angle -Omax? d'axe X2 // [320] (lignes de la figure 5b). Les positions de la 
deuxieme extremit6 sont simul6es dans le logiciel "jems" en d6sorientant le cristal de 
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r6f£rcnce par une rotation d'angle Gmax* d*axe X2 // [320] (lignes de la figure Sb). Unc 
rotation d'angle AO = 20niAx permet done de reproduire avec une bonne approximation, 
les largeurs AOg variables des bandes de Holz. 

Avant ramincissement par la technique TIB', le substrat de silicium 
est contraint, mais g6n6ralement peu d^formd. Apr&s amincissement, les contraintes 
engendrdes par Tinterface entre les parties a et b sont relachdes. En chaque point de 
coordonnde (yi, zj), les elements finis donnent done un ddplacement R(yi, Zi) de 
eomposantes u et v sur les axes yi et Z]. Ce ddplacement peut €tre ddcomposd en une 
translation loeale, une ddfomiation e(yi, zi)et une rotation locale d'axe X| et d*ang1e 
P(yb2)): 



1 


rav 


du 


~2 


1^1 





Les simulations dans le logiciel jams indiquent que ce sont ces 
rotations locales P(yi, Zj) le long du faisceau d'dlectrons (direction zo) qui sont la 
cause principale de I'dlargissement des lignes de Holz en bandes de Holz. Une rotation 
d'axe X| et d*angle P(yi, Z|) peut Sire ddcomposde en 3 rotations d'axe X2, d'angle 
0 (yi, Z|), d'axe y2, d'angle 0'(yi, Z|) et Z2, d'angle G'Xyi, Z|). Au premier ordre, seule la 
premifere de ces 3 rotations est importante et I'on a Tdgalitd : 



e(yi,z,)=0,98 p(yi,z,) 



En particulier, on a 6galement montrd que Tangle Omax est au premier 
ordre dgal a 0,98Pinax > ou pmax est la valeur maximale de Tangle p(y), zi) le long du 
faisceau d'61ectrons, c'est k dire le long de y^, ou ce qui revient au mSme a la valeur 
maximale de Tangle P(yi, Z|) le long de la direction de yi puisque notre lame est 
pdriodiqueenx) : 

Dans celte demonstration de faisabilit6, on utilise cette proprietd 
approchee pour optimiser les constantes du materiau. Chaque diagramme CBED pris a 
une profondeur zz est done caract6ris6 par un seul parametre A0(Z2) = 29max (za) 
(figure 3b). 

Les courbes de la figure 6 tracent les valeurs mesur6es 
exp6rimentalement et les valeurs 0,98* Ap., (avec Ap=2 *Pmax) d6terminees a partir 
des calculs par 616ments finis. 
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Les courbes calcul^es de la figure 6a et 6b uUlisent des param&tres du 
mal^riau NiSi pris dans la littdraturc. Pour reproduire au mieux les courbes 
exp6rimentales» on a fait varier les. paramitres du matdriau NiSi dans le calcul des 
courbes calculdes des figures 6c et 6d. Par essais/erreurs, on obtient les constantes de 
la couche NiSi (module de Young E et coefficient de Poisson v (en se plafant dans 
I'approxiniation isotrope). 

On donnera maintenant quelques pr6cisions sur rutilisation des 
61£ments finis. On a utilis6 un logiciel du commerce. On s*est plac6 dans 
Tapproximation des deformations planes. 

Trois 6tals particuliers sont importants. 

6tat 1 : On a suppos6 qu'il existe une temperature To=TiH-AT a 
laquelle le substrat de Silicium et la couche deNiSi seraient cohdrents, sans contrainte 
(cette temperature n'est pas forc6ment accessible, on pourrait la designer comme la 
temperature fictivc de coherence*). 

etat 2 : La lame mince d'epaisseur t, coherente a la temperature To est 
refroidie h la temperature experimentale Ti=To-AT Les contraintes sont en partie 
relaxees a la surface de la lame mince. 

On a utilise les lois classiques de la mecanique des solides, pour 
determiner les deplacements R(yi,2i). En particulier la relation entre contraintes et 
deformations est donnee par : 



0 + v)(l-2v) 
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ou AT = To-Tk E est le module d'Young, v le coefficient de Poisson, 
ct le coefficient de dilatation thermique. Ces relations s'appliquent dans le substrat de 
Si (Ea=1 56GPa, va=0,277 et aA==2,6e-6 K'^) et dans la couche mince de NiSi (Eb, Vb 
aB et AT ont ete optimises en minimisant la distance % ^^^^ '^s courbes 
experimentales et calculees (voir figure 6)) en utilisant Thypothese des deformations 
planes. Les deformations sont calcuiees dans le plan defini par la direction des 
electrons incidents yi // zO et Taxe 2|. 

Etat 3 : Avant extraction de la lame, le systfeme etait beaucoup plus 
contraint car les contraintes ne sont pas relSchees par les sxirfaces. Une fois les 
parametres (par exemple le coefficient d'eiasticite et le coefficient de dilatation 
tbermiques) evalues, on peut gr§ce au present proc6de, evaluer les contraintes dans le 
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circuit int6gr6 avant amincissement. Les conlrainlcs dans le transistor sont calculdes 
en utilisant les conslantes optimises : ce sont les contraintes prises au milieu d'une 
lame de tr6s grande ^paisseur toou celles d'une lame pdriodique d*^aisseur infinie. 
Remarques : 

5 Rl) Pour le point (iv), dans cette demonstration de faisabi1it6, on a 

ajuste manuellement les constantes d'dasticitd, mais une optimisation automatique est 
possible. 

R2) Pour le point (iv), on avons 6galement utilisd la relation 
10 approch6e : e,„ax = 0.98 pmax, c'est a dire que Ton n'a pris en compte que la largeur 
totale A9g des bandes de Holz. Un ajustage plus precis mais plus long consisterait ^ 
reproduirc les bandes de Holz dans leur ensemble. Chaque diagramme CBED ne serait 
plus alors caracl^ris^ par une seule valeur AG mais par un ensemble de profils 
identiques k ceux realises a la figure 7b. 
15 Ces profils repr^sentent en fait la fonction l-Ag^(s) ou Ag(s) est 

Tamplilude de la fonction d'onde diffract^e par les plans d'indice g et s mesure Tecart a 
la loi de Bragg pour ces plans g. Comme expliqu6 dans le livre pr6cit6 de Spence et 
Zuo, s et 9 sont relics par une formule simple et diff6rentes approximations peuvent 
etre r6alis6es pour calculer Ag(s)- On a montr^ que pour les lignes fines de Holz de 
20 diagrammes CBED pris dans la direction y2//[230], une approximation deux ondes ou 
une approximation cin6matique donnent le meme r^sultat physique. Dans 
Tapproximation cin^matique : 

A g (5, ) = £ exp(-2^ i gR( j^j , ^2 ) ^Vi -27risy^) dy^ 

25 

Le parametre a minimiser serait dans ce cas la somme des distances 
entre des courbes du type de la figure 7b et 7c. L4nt6ret de la formule pr6cedente est 
que non seulement les rotations, mais egalement les deformations sont prises en 
compte dans le terme R(y2,Z2) (m€me si on a montr6 que Teffet essentie] dans 
30 Telargissement des bandes de Holz est celui du aux rotations locales). 

R3) Les points (i) et (ii) ont 6t6 6galement realises sur r^chantillon 
n°2 d6crit dans la partie 3. 

R4) A difftrentes temperatures, le principe de la methode reste 
similaire et est relativement classique dans la mesure des deformations. Introduire des 
35 mesures a differentes temperatures revient a augmenter le nombre de donn^es 
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cxp^rimentales et 1e nombre de paramitres a optimiser : la m^ihodolo^e est tout h fait 
similaire. 
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REVENDICATIONS 
]. Proc6d6 de determination d'au moins une caract^ristique 
m^canique d'au moins une couche disposde k la surface d'un substrat caract^risd en ce 
quMI comporte : 

a) la realisation d'une lame (L) d'dpaisseur 1, notamment 
suffjsamment faible pour relficher les contraintes dans ladite 
lame et ayant deux faces sensiblement parallSles et dispos6es 
sensiblement pexpendiculairement i ladite surface de 
substrat ; 

b) la mesure sur ladite lame d*au moins un parametre de 
deformation du substrat k diff^rentes profondeurs par rapport 
k la surface ; 

c) la determination h partir au moins dudit paramfetre de 
deformation^ d*au moins une caracteristique mecanique de 
ladite couche. 

2. Procede selon la revendication 1, caracterise en ce qu'il comporte 
la realisation de plusieurs lames (Li, L2 -..) d'6paisseurs differentes ainsi que la mise 
en oeuvre de retape sur chacune desdites lames (Li, L2 ...)• 

3. Procede selon une des revmdications 1 ou 2, caracterise en ce que 
pour au moins une lame, retape est repetee a au moins deux temperatures 
differentes. 

4. Procede selon une des revendications precedentes, caracterise en 
ce que ladite mesure est realisee en generant pour des points du substrat situes a 
differentes profondeurs des diagrammes de diffraction d'un faisceau eiectronique 
convergent comportant des lignes de Holz. 

5. Procede selon la revradication 4, caracterise en ce que ladite 
determination comporte le releve de la largeur des lignes de Holz d'au moins certains 
desdits diagrammes pour au moins un plan cristallographique du substrat. 

6. Procede selon la revendication 5, caracterise en ce que Tetape c 
met en oeuvre a partir de la largeur desdites lignes de Holz correspondant audit plan 
cristallographique, le calcul pour chaque diagramme d*une rotation maximale p,nax 

le 

long de Paxe du faisceau d'6Iectrons. 

7. Procede selon la revendication 67, caracterise en ce que retape c 
met en oeuvre le trace d'au moins une courbe representant une dite rotation maximale 
en fonction de la profondeur h laquelle ont ete obtenus lesdits diagrammes. 
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REVENPICATIONS 

1. Proc&i6 de determination d*au moins une caract6ristLque 
m6canique d*au moins une couche disposSe h la surface d'un substrat caract6ris6 en ce 
qu*il comporte : 

5 a) la realisation d'une lame (L) d'6paisseur t, notanunent 

suffisanunent faible pour relScher les contraintes dans ladite 
lame et ayant deux faces sensiblement paralieies et dispos^es 
sensiblement peipendiculairement k ladite surface de 
substrat ; 

10 b) la mesure sur ladite lame d'au moins un parametre de 

defonnation du substrat k diff6rentes profondeurs par rapport 
k la surface ; 

c) la detemiination k partir au moins dudit parametre de 
deformation, d*au moins une caracteristique mecanique de 
15 ladite couche. 

2. Precede selon la revendication. 1, caracterise en ce qu'il comporte 
la realisation de plusieurs lames (Li, Li ...) d'epaisseurs differentes ainsi que la tnise 
en oeuvre de retapc b sur chacune desdites lames (Li, Li ...)• 

3. Precede selon une des revendications 1 ou 2, caracterise en ce jiue 
20 pour au moins une lame, I'etape b est repetee k au moins deux temperatures 

differentes. 

4. Precede selon une des revendications pr6c6dentes, caracterise en 
ce que ladite mesure est realisee en generant poxxr des points du substrat situes a 
differentes profondeurs des diagranunes de diffraction d'lm faisceau electronique 

25 convergent comportant des lignes de Holz. 

5. Precede selon la revendication 4, caracterise en ce que ladite 
determination comporte le releve de la largeur des lignes de Holz d'au moins certains 
desdits diagrammes pour au moins un plan cristallographique du substrat, 

6. Procdde selon la revendication 5, caracterise en ce que Tetape c 
30 met en oeuvre a partir de la largeur desdites lignes de Holz correspondant audit plan 

cristallographique, le calcul pour chaque diagranmie d'une rotation maximale Pmax. le 
long de Taxe du faisceau d'eiectrons. 

7. Precede selon la revendication 6, caracterise en ce que retape c 
met en oeuvre le trace d'au moins une courbe representant une dite rotation maximale 

35 en fonction de la profondeur k laquelle ont ete obtenus lesdits diagrammes. 
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8. Proc£d6 selon la revendication 7, caract6ris6 en cc que ]*itape c 
met dgalement en oeuvre par simulation le tracd de courbes reprdsentant la rotation 
maximale en fonction de la profondeur pour des valeurs possibles du module d*Young 
et/ou du coefficient de Poison, et la minimisation de Ticart entre les courbes simul6es 
el les courbes expirimentales pour determiner le module d'Young et/ou le coefficient 
de Poisson de ladite couche. 
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PROCEDE DE MESURE DE PARAMETRES PHYSIQUES D'AU MOINS UNE COUCHE D'UN MATERIAU AUX 
DIMENSIONS MICROMETRIQUES 


LE(S) DEMANDEUR(8) : 




COMMISSARIAT A L'ENERGIE ATOMIQUE 
31-33. rue de la Fdddratlon 
75015 PARIS 
FRANCE 


DESIGNE{NT) EN TANT QU'INVENTEUR(S) : 


D Norn 


ROUVIERE 


Pr6nonis 


Jean-Luc 


Adresse 


Rue 


1 1 , allSe de la Passto 




Code postal etville 


I3i8i2i4i0l MEYLAN (FRANCE) . 


Soci6t4 d'appartenance (fiicuUaHJ) 




U Norn 


CLEMENT 


Pr^noms 


Laurent 


Adresse 


Rue 


27, avenue Jeanne d'Arc 




Code postal et viile 


13 i8 i1 lO lO 1 GRENOBLE (FRANCE) 


Soci^te d'appartenance (facultatij) 




B Norn 


PANTEL 


Pr^noms 


Roland 


Adresse 


Rue 


2, rue Moucherotte 




Code postal et ville 


13 i8 13 12 lO 1 POISAT (FRANCE) 


Soci^td d'appartenance (factiUaHj) 




S'H y a plus de trois inventeurs, utilisez plusteurs formutafres. Indiquez en fiaut h droUe le de la page suhd du nomtire de pages. 


DATE ET SIGNATURE(S) 

DU (DES) DEMANDEUR(S) 

OU DU MANDATAIRE 

(Nom et qualit6 du slgnataire) ^ 




Paris, le 1 8 Juillet 2003 

Philippe JACQUARD - N** 92-4024 

Mandataire - CABINET ORES 





La lol n''78-17 du 6 Janvier 1978 relative k I'informatlque, aux fichteis et aux liberty s'applique aux r^ponses faites k ce formulaire. 
Elle garantit un droit d'accte et de rectification pour les donn^es vous concemant auprte de TINPI. 



PCT/FR2004/001877 




This Page is Inserted by IFW Indexing and Scanning 
Operations and is not part of the Official Record 

BEST AVAILABLE IMAGES 

Defective images within this document are accurate representations of the original 
documents submitted by the applicant. 

Defects in the images include but are not limited to the items checked: 

□ BLACK BORDERS 

□ IMAGE CUT OFF AT TOP, BOTTOM OR SIDES 
FADED TEXT OR DRAWING 
BLURRED OR ILLEGIBLE TEXT OR DRAWING 

□ SKEWED/SLANTED IMAGES 

□ COLOR OR BLACK AND WHITE PHOTOGRAPHS 

□ GRAY SCALE DOCUMENTS 

□ LINES OR MARKS ON ORIGINAL DOCUMENT 

□ REFERENCE(S) OR EXHIBIT(S) SUBMITTED ARE POOR QUALITY 

□ OTHER: 

IMAGES ARE BEST AVAILABLE COPY. 
As rescanning these documents will not correct the image 
problems checked, please do not report these problems to 
the IFW Image Problem MaUbox. 





